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Abstract 

We study static configurations of dark matter coupled to a scalar field responsible for the 
dark energy of the Universe. The dark matter is modelled as a Fermi gas within the Thomas- 
Fermi approximation. The mass of the dark matter particles is a function of the scalar field. 
We analyze the profile of the dark matter halos in galaxies. In this case our framework 
is equivalent to the model of the isothermal sphere. In the presence of a scalar field, the 
velocity of a massive object orbiting the galaxy is not of the order of the typical velocity of the 
dark matter particles, as in the conventional picture. Instead, it is reduced by a factor that 
quantifies the dependence of the dark matter mass on the scalar field. This has implications 
for dark matter searches. We derive new solutions of the Einstein equations which describe 
compact objects composed of dark matter. Depending on the scale of the dark matter mass, 
the size of these objects can vary between microscopic scales and cosmological distances. We 
determine the mass to radius relation and discuss the similarities with conventional neutron 
stars and exotic astrophysical objects. 

PACS numbers: 95.35. +d, 98.35. Gi, 98.80.Cq 
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1 Introduction 



Our present understanding of the history of the Universe assumes the presence of various contri- 
butions to its energy content. The total energy density is believed to take the critical value that 
results in a spatially flat Universe. The most accessible energy component is baryonic matter, 
which accounts for ~ 5% of the total energy density. A component that has not been directly 
observed is dark matter: a pressureless fluid that is responsible for the growth of cosmological 
pertrurbations through gravitational instability. Its contribution to the total energy density 
is estimated at ~ 25%. The dark matter is expected to become more numerous in extensive 
halos, that stretch up to 100-200 kpc from the center of galaxies [1]. The component with the 
biggest contribution to the energy density has an equation of state similar to that of a cosmo- 
logical constant. The ratio w = p/p is negative and close to —1. (For a review, see ref. [2].) 
This component is responsible for ~ 70% of the total energy density and induces the observed 
acceleration of the Universe [3] . 

The difficulty with explaining the very small value of the cosmological constant that could 
induce the present acceleration has motivated the suggestion that this energy component is time 
dependent [4, 5]. In the simplest realization, it is connected to a scalar field <j) with a very flat 
potential. The vacuum energy associated with this field is characterized as dark energy and 
drives the acceleration. If such a field affects the cosmological evolution today, its effective mass 
must be of the order of the Hubble scale, or smaller. 

It is conceivable that there is a coupling between dark matter and the field responsible for 
the dark energy [6]. In such a scenario it may be possible to resolve the coincidence problem, 
i.e. the reason behind the comparable present contributions from the dark matter and the dark 
energy to the total energy density. The construction of specific models is still in the early stages 
[7]-[ll]. Phenomenological difficulties arise because of the required flatness of the potential of 
the scalar field. The coupling between the field and the dark matter could lift the flatness of the 
potential through radiative corrections. We assume that at the level of the effective potential 
the required flatness is maintained, possibly guaranteed by some symmetry [12]. 

The presence of an interaction between dark matter and the scalar field responsible for the 
dark energy has consequences that are potentially observable. The cosmological implications 
depend on the form of the coupling, as well as on the potential of the field [7]. If the scale for 
the field mass is set by the present value of the Hubble parameter, then the field is effectively 
massless at distances of the order of the galactic scale. Its coupling to the dark matter particles 
results in a long range force that can affect the details of structure formation [13, 14, 15]. 

The attraction between dark matter particles mediated by the scalar field may lead to the 
formation of dense compact objects composed primarily of dark matter. We are interested in 
static solutions of the Einstein equations that describe such objects. Pure scalar field configu- 
rations, coupled to gravity, are known. An explicit analytical solution has been derived in ref. 
[16]. This solution, however, has a naked singularity at the origin. In the present work we show 
that the addition of fermions can eliminate the singularity and result in a configuration that has 
a natural physical interpretation. 

In the following section we develop the necessary formalism for our study. We assume that the 
interaction between dark matter and dark energy takes the form of an effective field-dependent 
mass term for the dark matter particles. We derive the Einstein equations and the equation 
of motion of the scalar field. In order to derive an equation of state for the dark matter, we 
model it as a Fermi gas in the Thomas-Fermi approximation. In section 3 we summarize known 
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solutions of the system of equations of section 2. We also discuss a class of solutions that gives 
an approximate description of galaxy halos in the presence of the scalar field. In section 4 we 
derive new solutions that describe compact objects composed of dark matter and held together 
by the scalar interaction. Section 5 contains a summary and our conclusions. 



2 Formalism 

2.1 The basic expressions 

We assume that the dark matter consists of a gas of weakly interacting particles. The mass m 
of the particles depends on the value of a slowly varying classical scalar field <f> [6]. For classical 
particles, the action of the system can be written as (see ref. [8] and references therein) 

S = j d'x^-g [m 2 R - \^^!T ~ U{4>)) - E / (2-1) 



with dri = y — g^ v (xi)dx^dx^ and the second integral taken over particle trajectories. Variation 
of the action with respect to 4> results in the equation of motion [8, 17] 



y/—g dx^ V dx 
where 

W = ^=Y,jdn (x - Xi ) = - dln y )} T*. (2.3) 

The energy-momentum tensor associated with the gas of particles is 

dxf dx v 

^ > / (It, ///(o(.r,)) 



T^ = -LE/ *i ^-^5^\x-x t ). (2.4) 

J—g J dn dn 



In the following we shall look for solutions of eqs. (2.2), (2.3) employing an approximation for 
the form of the energy-momentum tensor T^ v . We shall assume that it takes the diagonal form 
T% = diag(-p,p,p,p). 

We are interested in stationary, spherically symmetric configurations of the system. We 
consider a metric of the form 

ds 2 = -B(r)dt 2 + r 2 {d8 2 + sin 2 6 # 2 ) + A{r)dr 2 . (2.5) 

The Einstein equations are 

11 1 I A' 1 | 

T 2 ^ ~ T 2 ~ rA 2 ~ 2M 2 \ 2A 



I A' IB' A'B' B' 2 B" 



~2rA 2 + 2r AB AA 2 B AAB 2 ' 2AB 



- U(0) - p) (2.6) 

+ ^b = ^(-hf-™**) (2 - 7) 



11 1 IB' 
^ 2 A~^ 2 + rB~A 2M 2 
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where a prime denotes a derivative with respect to r. The scale M is defined as M = (167rG) 1 I 2 , 
where G is Newton's constant. The equation of motion (2.2) for the field <f> becomes 



'2 A' B' 
+ ' r ~ 2A + 2B 



b' = A 



dU din-mid)) , „ ' 



(2.9) 



Combining eqs. (2.6), (2.8), (2.9), we can eliminate A' and B' from the equation of motion of 
the field 6 



l + A 



A 



dU dlnm( 



(p-3p) 



(2.10) 



It is more convenient to replace one of the equations of motion by the conservation of the energy- 
momentum tensor. Combining the only non-trivial conservation equation with the equation of 
motion (2.9) gives 

dp ^ x d\nm(4>) dd> , % 1 dB 



dr = ^- p) ^f L Tr- {p + P) 



2B dr 



(2.11) 



We also give some other useful forms of the Einstein equations. Eq. (2.6) can be integrated 
with the result 



where 



A(r) 

dM(r) 
dr 



1 



1 M{r) 



8vrM 2 



Eqs. (2.6), (2.8) can be combined in order to obtain 

dln(BA) 1 



k'2 



+ A(p + p) 



r. 



dr 2M 2 

The Newtonian limit is taken easily if we combine eqs. (2.6)~(2.8) in order to obtain 

1 



(2.12) 
(2.13) 

(2.14) 



B" + -B' = -r^AB 
r 4M 2 



(1 + A)p + (5 + A)p -4U + 



A 



-10 



+ 



8M 4 



ABr A 



1 



h'2 



A {p - U){ P + p) + cf> u \-{p + 3p) + ^ - U 



(2.15) 



For B = 1 + 2$, with |$| <C 1, we can set A = B = 1 in the r.h.s. of the above equation and 
neglect the second line. In this way we find 



AM 2 



(p + 3p- 2U) 



(2.16) 



for the Newtonian potential. In the same limit, eq. (2.10) becomes independent of the gravita- 
tional field and reads 

„ 2 , dU d\a.m{4>) 



^ + ^ = ^ + =^(p-3p). 



(2.17) 
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2.2 The Thomas- Fermi approximation 



The most popular models of dark matter assume that it consists of a gas of weakly interacting 
fermions. In the previous section we derived the formalism for our study neglecting the quantum 
nature of the particles. As we are interested in dense dark matter configurations, we must 
account for the fact that the quantum properties may become relevant at sufficiently high density. 
In our treatment we shall only take into account the physical implications of the exclusion 
principle. We neglect particle scattering, but we assume that the gas is described by a Fermi- 
Dirac distribution. This is the essense of the Thomas-Fermi approximation. 

In a local frame at every point in space the particles are described through the distribution 



f(p) 



cxp 



^/p 2 + m 2 ((j)(r)) — u(r) 



+ 1 



(2.18) 



The chemical potential fx and temperature T, as measured by a local observer, are functions of 
the radial coordinate r, while the mass of the particles depends on the local value of the field (f>. 
The pressure, number density and energy density are given by 



V 



cxp -- 



\Jp 2 + m' 2 (4>(r)) — u(r) 



T(r) 



+ 1 



n = 



1 

471-3 
1 

471-3 



d 3 p f(p) 
d 3 p f(p) 



3\/p 2 + m 2 {4>{r)) 



= ^3 J ' d 3 p f(p) \Jp 2 + m 2 



(r)). 



The above quantities are related through 

p = fin — p + Ts, 

where s is the entropy density. An important identity is 

dp d\nm(4>) d\nm{ 
~ — 71 O ~ P) = 



rpfl 



d(f) d(f> ' ' ' d(j) 

Combining this identity and the conservation equation (2.11) results in 

dp . . 1 dB 



(2.19) 

(2.20) 
(2.21) 

(2.22) 
(2.23) 

(2.24) 



Evaluating the r-derivative of p, given by eq. (2.19), and employing the relation (2.22) permits 
us to rewrite eq. (2.24) as 



dr 

This equation is satisfied if 

T(r) = n/jBir), 



dr 



K r ) = Po/\/B(r). 



(2.25) 



(2.26) 
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These are the standard expressions that describe the behaviour of a Thomas-Fermi fluid in a 
gravitational field. 

The derivation of the equation of motion (2.9), on which our discussion is based, neglected 
the quantum nature of the particles. However, within the Thomas-Fermi approximation the 
only modification relative to a gas of classical particles is the change of the distribution. When 
this is taken into account in the energy density and pressure, eq. (2.9) remains valid. For T = 0, 
this has been shown explicitly in ref. [18], starting from the equation of motion of a fermionic 
field (the Dirac equation) in a curved background. 

A consistent solution can be obtained by integrating eqs. (2.6), (2.8), (2.9) or (2.10), (2.11). 
Eq. (2.7) is then automatically satisfied. As we have shown, eq. (2.11) can be replaced by eqs. 
(2.26). In this way, the system of equations to be solved is reduced to eqs. (2.6), (2.8), (2.9) or 
(2.10). 

3 Solutions 

3.1 Compact astrophysical objects 

The Thomas-Fermi approximation in a gravitational background has been considered repeatedly 
in the past. In refs. [18, 19] the formalism was developed for a degenerate fermionic gas with 
T = 0, interacting with a scalar field through the mass of the particles. In ref. [20] a fermionic 
gas with non-zero T and [i, but no scalar field, was considered. In the previous section we 
developed the formalism for the more general case of a fermionic gas with non-zero T and //, 
interacting with a scalar field through the mass of the particles. 

The set of equations we derived has various solutions that depend on the assumptions made 
about the matter content and the boundary conditions. We summarize briefly known solutions 
of eqs. (2.6), (2.8), (2.9) or (2.10) with T, p obeying eqs. (2.26). 

a) In the absence of a scalar field, a degenerate fermionic gas at zero temperature can form a 
compact astrophysical object, in which the gravitational attraction is balanced by the degeneracy 
pressure. Such solutions have been discussed in ref. [20]. The binding energy that keeps the star 
together is provided by gravity. For this reason these solutions are not stable for small values of 
the total fermionic number. The mass is a decreasing function of their radius for small masses. 
The possibility of non-zero temperature T <C \i has also been considered. 

b) If a scalar field with an appropriate potential is considered, a degenerate fermionic gas can 
form astrophysical objects that have been termed fermion soliton stars [18] or fermion Q-stars 
[19]. These objects have T = and the fermions can be assumed to be ordinary nuclei. In 
particular examples, the scalar field has been identified either with the cr-field associated with 
the nuclear forces, or with a mesonic condensate (ir- or i^-condensate) . The main difference with 
the previous class of solutions is that these objects are bound even in the absence of gravity. The 
value of the scalar field that minimizes the fermionic mass does not coincide with the mimimum of 
the potential. A new phase exists in which the fermions are essentially massless, while the scalar 
field is displaced from its vacuum expectation value. An appropriate choice of the couplings 
guarantees that this phase has positive binding energy. The gravitational interaction increases 
the pressure at the center of the object, as well as the binding energy per fermion. 

The method for constructing solutions of this type has been described in ref. [19]. In the case 
that has been considered more often, the fermionic mass is m{4>) = g<fi- The chemical potential 
fj, is taken in a range such that the energy per fermion in the interior is smaller than the mass 
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of a free fermion. As a result, the fermions are trapped in the interior. A static, spherically 
symmetric solution of arbitrarily large radius can be constructed by making use of eq. (2.23) in 
order to interpret eq. (2.17) as the equation governing the motion of a particle in the potential 
Ueff = p — U. This is possible through the replacement r — ► t and (f> — > r. The theory is assumed 
to have a vacuum characterized by a field value (f> = 4>i for which {7 = and p = p = 0. If the 
gravitational field is neglected, a new solution <j> = </>q is needed, for which d(p — U)/d(f> = and 
p — U = 0. This can be achieved by adjusting the chemical potential. If d(p — U) / dcj) = — e at 
0o with e — ► 0+, any solution with (j>'{r = 0) = takes a long "time" r before moving away from 
00 towards larger values. At r ~ R the field starts "rolling" quickly towards 4>\. The "friction" 
term ~ 2/r has a negligible effect, so that the "rolling" takes place between the two maxima 
of U e ff near 4>q and <f>i without loss of energy. Soon after r = R the fermionic density becomes 
zero, as m{4>) exceeds p. The interior of the configuration corresponds to r ^ R. 

After the inclusion of gravity, the solution in the interior is again determined by the condition 
d(p — U)/dcf> = [19]. The r-dependence of p, according to eq. (2.26), implies that the field 
4> now varies with r as well. The variation is very slow, so that the l.h.s. of eq. (2.10) can be 
approximated by zero. The quantity p — U is non-zero in the interior of the solution [19]. Its 
value at the center must be adjusted so that the exterior solution, which involves only the scalar 
field, is matched in the presence of gravity. The mass of these soliton stars is an increasing 
function of the radius for small masses. 

c) Pure scalar field configurations, coupled to gravity, have been discussed in ref. [16], where 
an explicit analytical solution has been derived. This solution, however, has a naked singularity 
at the origin. The addition of fermions is necessary in order to eliminate the singularity. In 
section 4 we present a solution in which this is achieved. 

3.2 Dark matter halos 

The most prominent feature of the distribution of dark matter particles in galaxy halos can be 
understood using the formalism of section 2. The dark matter distribution results in approxi- 
mately flat rotation curves 1 for objects orbiting the galaxies at distances r ^ 10 kpc [1, 21]. An 
analytical understanding of this behaviour is possible within simple models of the dark matter 
gas, such as the isothermal sphere. If these simple models are extended through the addition of 
a scalar field to the theory, an analytical treatment is still feasible. For our study we model the 
dark matter gas as a Thomas- Fermi fluid, described by the formalism of section 2. We assume 
that the gas is non-relativistic and non-degenerate, with p <C p. In the absence of a scalar field 
{4> = 0, U = 0, m = mo) we can define the non-relativistic chemical potential as po = po — mo- 

x It must be emphasized that the rotation curves deduced from observations or computer simulations are not 
exactly flat, but reveal a structure related to the underlying physics of structure formation [22]. In particular, the 
distribution of matter in the galaxy cores is an open problem, for which analytical approaches are not currently 
available. Theoretical studies rely mainly on computer simulations. Here we make the assumption that the 
approximate flatness of the rotation curves outside the galaxy cores is related to the virialization of the dark 
matter gas. The simple model of the isothermal sphere captures the essence of this assumption, while it can 
be easily extended to include the scalar field. The solution of the Einstein equations that we summarize is 
equivalent to the model of the isothermal sphere. It cannot explain the details of the observed or simulated 
rotation curves, apart from the gross feature of the approximate flatness outside the cores. However, it leads to 
important conclusions for the dispersion of the velocity of the dark matter particles, which are crucial for dark 
matter searches. 
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The number density is 



n ~ 2 (^) 3/ expf-^ + ^ ~n exp (-«$), (3.1) 



V 2vr J \ T T 

with 

/m r \ 3/2 fp \ m . . 

The validity of the above expressions requires /2o < and Tq <C mo, T <C | A*o I - The gravitational 
field can be studied in the Newtonian limit. Then, eq. (2.16) can be written as 

d 2 u 2du 

^ + ^ + eXpU = °' (3 ' 3) 

with u = —a§, z = (3r, (3 2 = apo/AM 2 , and po — mo n o- The solutions of this equation that 
are regular at z = behave as expu = 2/z 2 for large z [1]. For large r the Newtonian potential 
varies only logarithmically with r, while the integrated mass of the dark matter scales linearly 
with r. This leads to flat rotation curves for objects orbiting the galaxy [1]. 

In ref. [23] the possibility was considered that the mass of the dark matter particles depends 
on the field (f> that varies slowly with the radial distance r. The analysis can be easily reproduced 
through the formalism we developed in section 2. We approximate the number density of the 
dark matter particles by the expression (3.1), with mo replaced by m(</>). As before, we work 
in the weak-field limit, in which B(r) ~ 1 + 2$(r), with |$(r)| <C 1. We assume that the 
field 4> is displaced from its asymptotic (r — > oo) value only by a small amount, so that the 
approximation m(4>) = m(0o) + [dm((j>o) / dcf)] 5<p = mo + m' 5(f> can be employed. We identify 
0o as the value of the field at the center of the galaxy (r = 0). We work in the leading order 
in 5(f), and assume that m'/m ~ m /mo for all r. Our treatment is relevant up to a distance 
r\ ~ 100-200 kpc beyond which the dark matter becomes very dilute. For r ^ ri, we expect 
that (p quickly becomes constant with a value close to <f)(r\) = 4>\. This is the value that drives 
the present cosmological expansion [4, 5, 6]. Here we assume that the cosmic evolution of (f>± is 
negligible for the time-scales of interest, so that the asymptotic configuration is static to a good 
approximation. 

Within the leading order in 5(f>, we can approximate dU/d<p by a constant between r = 
and r = oo. For the scalar field to provide a resolution of the coincidence problem, the two 
terms in the r.h.s. of eq. (2.2) must be of similar magnitude in the cosmological solution. This 
means that dll/d(p must be comparable to (m /mo)p 00 . We expect to be a fraction of the 
critical density, i.e. p^ ~ 3 keV/cm 3 . On the other hand, for the spherically symmetric solution 
the energy density in the r.h.s. of eq. (2.2) is that of the galaxy halo (~ 0.4 GeV/cm 3 for our 
neighborhood of the Milky Way). This makes dU/d<p negligible in the r.h.s. of eq. (2.2) for a 
static configuration. The potential is expected to become important only for r — > oo, where the 
static solution must be replaced by the cosmological one. Similar arguments indicate that we 
can neglect U relative to p. The scale for the field mass is expected to be set by the present value 
of the Hubble parameter. Then the field is effectively massless at distances of the order of the 
galactic scale. The above indicate that, if the deviation of the scalar field from its cosmological 
value is small, the form of the potential plays a negligible role at the galactic level. For this 
reason we set U = in our discussion. 
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The scalar field generates a new long-range scalar interaction, whose strength relative to the 
gravitational interaction is determined by the parameter 



K 2 = 4M 2 (^°) . (3.4) 

If the new interaction is universal for ordinary and dark matter, the experimental constraints 
impose k <C 1. In this case, it is reasonable to expect a negligible effect in the distribution of 
matter in galaxy halos. However, if <p interacts only with dark matter this constraint can be 
relaxed. In the following we assume that the scalar field interacts only with the dark matter. 

As before, we assume p <C p and employ the nonrelativistic chemical potential po = po — mo- 
The number density of dark matter can be written as 

n ~ 2 \— g^- J exp y ^ + £ J ~ n exp (-a* - aS<f>) , (3.5) 

with no, a given by eqs. (3.2) and 

<ty = 0-fo = 0-#r = O), « = ^- (3.6) 

Jo 

The energy density of dark matter at the center of the galaxy is po = niQUQ. 

It is important to emphasize that the assumption that the number density is given by eq. 
(3.1) does not require the presence of thermal equilibrium. In an alternative approach, followed 
in ref. [23], the dark matter can be considered as a dilute, weakly interacting gas with an energy- 
momentum tensor T^ u = diag(— p, p,p,p). Motivated by the model of the isothermal sphere [1], 
we can assume that p(r) = p(r) {v 2 d ) = m((f>(r)) n(r) (v^), with a constant velocity dispersion 
(v 2 ,) < 1. In the weak field limit and for p <C p, the conservation of the energy-momentum 
tensor (2.11) gives 

j/ = -p#-p<(6<l>y. (3.7) 
m 



Integration of this equation results in 

n ~ no exp I — —^r — —j-jr ) ■ (3.8) 



$ 777-0 ^ 



By defining an effective temperature To through the relation (v 2 ,) = To /mo, we reproduce eq. 
(3.5). This indicates that the assumption of thermal equilibrium is not required for the emergence 
of eq. (3.5). The parameter To that appears in the various expressions of section 2 does not 
correspond necessarily to the physical temperature. In many cases of physical interest it is 
simply a measure of the typical velocity of the Fermi gas. 

In the non-relativistic, weak-field limit, with p = U = 0, eq. (2.16) becomes 

$,/ + - $ ' = 7TJ2/ 9 oexp(-a$-a^). (3.9) 

Similarly , eq. (2.17) becomes 

{5(t))" + -m)' = ^p exp (-a* - &S</>) . (3.10) 
r mo 
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A linear combination of eqs. (3.9), (3.10) results in eq. (3.3) where now 



u = -a$ - aScj), z = fir, (3.11) 

with 

o2 ( _a . ~m' \ 2 m p 

= [up + a ^r ) p0 = (1 + K ] n ii^- (3 - 12) 

The solutions that are regular at z = approach the form 



u = In 



2 \ 1 



z 2 1 s/z 



d\ cos I m z I + «2 sm I In z 



+ ... (3.13) 



for large z. Another linear combination of eqs. (3.9), (3.10) gives 



d 2 v 2dv , 



with 



,m' n 



v = a5(j)-AM 2 ^-a^ = -K 2 a^ + a5(t). (3.15) 

The solution of this equation is v = cq + c\j z. 

The velocity v of a massive baryonic object in orbit around the galaxy, at a distance r from 
its center, can be expressed as 

r3>' z ( du dv\ , 



where 



The asymptotic form of n(z), ■u(z) indicates that v ~ v c for large z. The dominant correction 
to the leading behaviour arises from the term ~ ^/\fz in eq. (3.13). The function v(z) gives a 
higher order correction. This means that the presence of the field eft is not expected to cause 
significant modifications to the shape of the rotation curves relative to the = case. For 
a vanishing field the rotation curves are again governed by the solution of eq. (3.3), given by 
eq. (3.13). This simple analysis indicates that the approximately flat rotation curves outside 
the galaxy cores are a persistent feature even if the dark matter is coupled to a scalar field 
through its mass. However, numerical simulations are probably necessary in order to reproduce 
the detailed form of the curves. 

We can use eq. (3.17) in order to fix To/too = (v^) for a given value of k. The effect of the 
scalar field is encoded in the factor k 2 = 4M 2 (to /too) 2 . When this is small, the velocity of an 
object orbiting the galaxy is of the order of the square root of the dispersion of the velocity of 
the dark matter particles. If n 2 is large the orbital velocity can become much smaller than the 
typical dark matter velocity. This behaviour persists even if baryonic matter is added near the 
center of the galaxy [23]. Within our extension of the model of the isothermal sphere, the flatness 
of the rotation curves is largely insensitive to the details of the matter distribution within the 
galaxy core. 

The scenario we considered may be more interesting for dark matter searches than the 
conventional one. For large values of k the velocity of dark matter particles exceeds significantly 
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the observed rotation velocity (~ 220 km/s for the Milky Way). The estimated local energy 
density of dark matter remains the same as in the case with re 2 = 0. It is ~ 0.4 GeV/cm 3 for 
our neighborhood of the Milky Way. As a result the flux of dark matter particles towards a 
terrestrial detector is larger roughly by a factor (1 + re 2 ) 1 / 2 relative to the re 2 = case. 

A detailed calculation of the counting rates in detectors must take into account the motion 
of the Earth around the Sun and the motion of the Sun through our galaxy [24]. The velocity 
distribution of dark matter is 

/(v)~exp(-!^±|£±MJ, (3 . 18) 

where the magnitude of the velocity of the Earth relative to the Sun is |ve| ~ 30 km/s, and 
that of the Sun relative to the galactic rest frame |v*s| ~ v c ~ 220 km/s. For re ^ 1, we have 

(v 2 ,) 3> |ve|, |vg|. The motion of the Earth and the Sun are expected to give only a small 
modification of the dark matter flux towards the Earth. As a result, the seasonal variation of 
a possible dark matter signal decreases for increasing re 2 . Typically, the cross section for the 
elastic scattering of halo particles by target nuclei is independent of the particle velocity for very 
low velocities [26, 21]. The leading effect of a non-zero value of re is that the counting rates, 
that are proportional to the dark matter velocity, are increased by the factor (1 + re 2 ) 1//2 . This 
makes the dark matter easier to detect. Existing bounds on dark matter properties from direct 
searches can be extended to include the case of non-zero re. The bound on the cross section for 
the interaction of dark matter with the material of the detector must be strengthened by the 
factor (1 + re 2 ) 1 / 2 . 

Corrections to the above simple picture may arise for large values of re 2 . For example, the 
nuclear form factors that must be included in the calculation of the counting rates have a velocity 
dependence [21]. If the velocities are larger than the conventionally assumed ^ (v 2 ,) ~ 220/^/2 
km/s, these form factors may result in significant modifications of the expected rates for the 
various experiments. A separate detailed study is necessary in order to address this issue. 

The allowed range of re is limited by the observable implications of the model that describes 
the dark sector. The dependence of the mass of dark matter particles on an evolving scalar field 
during the cosmological evolution since the decoupling is reflected in the microwave background. 
The magnitude of the effect is strongly model dependent. In the models of ref. [7, 9] the 
observations result in the constraint re 2 ^ 0.01. In the model of refs. [8, 15] the scalar interaction 
among dark matter particles is screened by an additional relativistic dark matter species. As a 
result, the model is viable even for couplings re 2 ~ 1. A similar mechanism is employed in ref. 
[10]. In this model the interaction between dark matter and dark energy becomes important 
only during the recent evolution of the Universe. In general, an interaction that is effective for 
redshifts z 1 — 2 is not strongly constrained by the observations. 

Before concluding this section we would like to give a more physical interpretation of the 
main physical effect. The dark matter in galaxy halos is expected to be a virialized gas, whose 
velocity is determined by the potential. The presence of the field associated with the dark 
energy adds an attractive scalar potential to the gravitational one. Their relative strength is 
given by the parameter re 2 defined in eq. (3.4). As the scalar field is essentially massless, both 
potentials have the same radial dependence. Their combined effect results in a total potential 
that is stronger by a factor 1 + re 2 compared to the conventional case. As a result, the kinetic 
energy of the virialized gas and the velocity dispersion are inreased by the same factor. 
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4 New astrophysical objects 



4.1 The ingredients for a solution 

In the previous section we discussed the effect of the scalar field on the distribution of dark 
matter in galaxy halos. Our basic assumption was that the field is not displaced significantly 
from its value at r — > oo. In this section we are interested in solutions with larger deviations 
of the field from its asymptotic value. For this reason we need a treatment that goes beyond 
linearized gravity. 

In the framework we are considering, the potential U((f>) does not have a minimum. The 
typical case, which we shall take as a working example, is an exponential potential U (</>) = 
Uo exp(— c<f>/M) [4, 6, 7]. There are two possibilities, depending on the form of m(<f>): 

a) If the mass m is a decreasing function of <j) (such as the decreasing exponential of refs. [6, 7]), 
the pressure p becomes an increasing function of (f>. This is obvious from eq. (2.23), taking 
into account that p < p/S. For a weak gravitational field we can make use of eq. (2.23) in 
order to interpret eq. (2.17) as the equation governing the motion of a particle in the potential 
U e ff = p — U. This is possible through the replacement r — > t and <j> — > r. The potential U e ff 
is an increasing function of (p. Any static configuration must start with a large value of (j) at 
r = 0, which then "rolls" for r — > oo to a smaller value <\>\ that is relevant for the cosmological 
solution. In more physical terms, this behaviour can be understood through the tendency of 
the system to minimize the energy by reducing the particle mass. As this mass vanishes for 
4> — ► oo, we expect that the field becomes larger than its cosmological value in the interior of the 
most stable static configuration. If we assume that during the cosmological evolution (which we 
neglect in the study of static configurations) the field (j) moves from smaller to larger values, such 
configurations may emerge dynamically if the field in certain regions grows faster than in the 
bulk of space. A strong gravitational field is not expected to induce qualitative modifications of 
this behaviour. 

b) There are also scenaria in which the dark matter mass is an increasing function of 4>. For 
example, in refs. [8, 15] the mass of the fermions that compose the dark matter is determined 
by the expectation value of the scalar field through a Yukawa coupling. We do not consider all 
the details of the scenario of refs. [8, 15], which involves two dark matter families with different 
masses. Instead, we limit ourselves to one family of dark matter particles. If the mass is an 
increasing function of the field, it is reasonable to imagine that, during the evolution of the field 
from smaller to larger values, certain regions of space become disconnected from the general 
evolution because dark matter gets trapped there. The field remains close to zero in these 
regions. The vanishing of the dark matter mass for zero <p makes such a scenario energetically 
favourable. 

In the examples we shall consider we shall assume that the dark matter mass is linear in <fi. 
The essence of our assumption is that the mass vanishes for some finite value of (f>, which can be 
set to zero by an appopriate field shift. The potential U{4>) retains its form, as the field shift can 
be absorbed in the pre-exponential factor. The strong dependence of the mass on 4> causes the 
field near dense concentrations of dark matter to deviate strongly from its asymptotic value 
The field tends to approach zero in the center of such concentrations, so that the dark matter 
particles become massless there. Our framework corresponds to the second case discussed above. 
However, the same physical mechanism could also operate in the first case, for which the field 
is expected to be larger than (f>\ in the interior of concentrations of dark matter. Astrophysical 
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configurations similar to the ones we describe in the following are expected to appear in this 
case as well. 

The solutions we shall derive are very different from the solutions of the previous section, 
as the dark matter gas will be degenerate (T = 0). It is possible to find solutions with T / 
following an analogous procedure, even though the technical difficulty may be greater. 

We have seen that the effective temperature of the fermionic gas can be related either to the 
velocity dispersion of a weakly interacting collection of particles (as for galaxy halos), or to the 
real temperature of a thermalized gas (as for the interior of stars). Our formalism can describe 
static solutions in both cases. The value of the temperature is determined by the ability of a 
collapsing gas to lose momentum, and the possible interactions between the particles that can 
establish thermalization. We shall not address this dynamical issue here. In the remaining of 
this section we shall present static solutions that describe a system very different from the galaxy 
halos we discussed in the previous section: a completely degenerate fermionic gas, interacting 
with a scalar field. 

The pure scalar field configuration that solves the Einstein equations in the absence of a 
potential and dark matter has been discussed in ref. [16]. We expect this solution to be realized 
for distances r > R, for which the dark matter density becomes negligible (of the order of the 
critical density of the Universe). The dark matter gas is concentrated in the region r < R. The 
radius of the compact object is R. For U = and large r, the leading terms of the pure field 
solution read 

<f> = fa - >yM^ B = A~ 1 = l-^±. (4.1) 

r r 

The free parameter R s determines the mass of the object as seen from an observer at infinity 

M tot = 8itM 2 R s . (4.2) 

The asymptotic value of the field for r — > oo is fa, while the free parameter 7 determines the 
derivative energy density 

--^-^ <«> 

As long as K, 3> U, the potential can be neglected. For distances r of the order of the Hubble 
radius this condition is not expected to be satisfied. At such scales the solution we are describing 
must be replaced by the time-dependent cosmological solution. 



4.2 The equations 

We would like to modify the solution of ref. [16] by adding dark matter in the region r < R. 
This addition can eliminate naturally the singularity that appears in the pure field configuration. 
The simplest form of the mass of the dark matter particles is 

m(<f>) = m + g<t>, (4.4) 

where we have assumed a Yukawa interaction between the dark matter fermions and the scalar 
field. The mass mo is expected to be associated with a new scale a, different from the Planck 
scale M. For example, a could be set by the expectation value of a heavy Higgs field in a Grant 
Unified Theory. In order for the Yukawa interaction not to generate a contribution to the mass 
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much larger than tjiq, the Yukawa coupling must be taken g S5 a/M. Then the constant tuq can 
be absorbed in a shift of <\> of order M. In the following we shall assume that the mass is 



m(<f>) = a^, (4.5) 

where we have taken g = a/M for simplicity. It must be pointed out that our assumption does 
not provide an explanation for the smallness of the Yukawa coupling. The appearance of a small 
parameter is inevitable in a system with two physical scales: the expectation value of </>, taken 
of order M, and the physical mass of the dark matter particles, often assumed to be as low as 
the TeV range. The situation is similar to the flavour problem in the Standard Model, for which 
a deeper understanding is lacking. 

We assume that the scalar field has a potential of the form 

U(<P)=Ca*e W (-c^y (4.6) 

with c = 0(1). The value of 4> that is relevant for the present cosmological evolution is given 
by the requirement that U((f>) be of the order of the critical energy density U (0) ~ 10~ 47 GeV 4 . 
This value 



M V c 



108 + lnC + 41n 



a 



(4.7) 



.GeV, 

must be approached by our static solution for r — > oo. A possible constant contribution m* to the 
mass defined in eq. (4.5) has been absorbed in <f>. Through this redefinition <f>/M — ► 4>/M+m*/a. 
We assume implicitly that = 0(c) in order to avoid unnaturally large values of 0. 

The equations of motion become more transparent if we define the dimensionless variables 

All other dimensionful quantities are multiplied with appropriate powers of a only, in order to 
form dimensionless quantities denoted as tilded. We also define the quantity 

6 = I = (4 - 9 > 

Mo po 

We have the relations 

m{4>) = 4>, £(f) = -fi=. (4.10) 

y/B(f) 

We define the surface of the compact object as the point at which the fermionic density and 
pressure vanish. If the surface corresponds to a value f = R = 0(1), the physical radius is given 
by the relation 

"sW^V 1 *. (411) 

km \ a J 

Similarly, from eq. (4.2), we find 

M tot /GeV\" 



Mr: 







-0.115^— j R s . (4.12) 
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Another important characteristic of the solution is the total fermionic number, which we assume 
to be conserved. The fermionic number density j° is the time component of a covariant 4- vector. 
We can deduce its value from the local density n, through the appropriate tetrad factor V a , 
where g^ u = V^V^r]^. The total fermionic number is 



N 



ri gd 6 xf = I 47rr 2 drVAn 
Jo 



MY 
o J Jo 



4-irf 2 dfyfAh = { — J N. 



The equations of motion become 
~1 + A 1 



2 Af Uj(4>) + -(p-p) 



= A 



dU 1 , 

-j + T (p - 3p) 



11 1 1 A' 1 / 1 , .. 



f A 2 



2A 



11 1 1 B 1 1/1 
where a prime denotes a derivative with respect to r. We also have 



n 



P 



1 /_ 2 ~2\ 3 / 2 

-m j 



/i + \/li 2 — fh 2 



(4.13) 

(4.14) 
(4.15) 
(4.16) 

(4.17) 
(4.18) 

(4.19) 
(4.20) 

We need four initial conditions for the system of equations (4.14)-(4.16). One of them is 
imposed by the regularity of the spherically symmetric solution at f = 0: <fi'(0) = 0. Another 
one is implied by eqs. (2.12), (2.13): A(0) = 1. The value of B(0) can be chosen arbitrarily. 
However, the normalization of eqs. (4.1) implies that AB(r — > oo) = 1. This means that 
AB(f — ► oo) = (po/a)" 2 , where we have used the definition (4.9). As a result, the choice 
of B(0) determines the chemical potential. Finally, <^>(0) must be chosen so that <j)(f — ► oo) 
reproduces correctly the present value 4>\ of the scalar field in the cosmological solution, as 
given by eq. (4.7). (We assume that the time scale of the cosmological solution is very large 
and neglect the time dependence of <p(r — > oo).) 



3vr 2 
1 

24vr 2 
1 

8^ 



h 2 ( 2p 2 - 5m 2 ) + 3m 4 In 

~ nr ~^,( n ~2 ~2\ ~4, ( p + ^Jp 2 -m 2 

H\j fi — m z ylfi — m J — m In I 



for jl > m, and h = p = p = for p, < fh. Finally, 

U{4>) = Cexp(-c^) . 



4.3 The solution 

In fig. 1 we display the form of the solution for a model with cr/GeV = 1. The chemical potential 
takes the value Ho/a ~ 0.33, while the scalar field approaches the value 4>i/M = 1 for large r. 
The potential U(4>) has c = 1, C ~ 3 x 10~ 47 and is negligible in the range of distances of 
interest. We plot the quantities 4>, p, A, B as functions of f. We observe that the scalar field 
approaches zero near the center of the solution, so that the fermions become almost massless 
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Figure 1: The solution of eqs. (4-14)> (4-15), (4-16) for a model with a/GeV = 1, Ho/o~ — 0.33, 
<j)\/M = 1, c = 1, C ~ 3.4 x 10~ 47 . We plot (f>, p, A, B, as a function of f. The object has 
R ~ 46, R s ~ 17, N ~ 920. 



there. The pressure and density of the fermionic gas vanish for f > R ~ 46. The value of R 
determines the radius of the astrophysical object through eq. (4.11): R ~ 16 km. The metric 
components A and B deviate significantly from 1 in the interior of the solution. For this reason a 
fully non-linear treatment of the Einstein equations has been necessary. The mass of the object 
can be deduced from the asymptotic form of A or B for f — > oo through the second of eqs. (4.1). 
We find R s ~ 17, which corresponds through eq. (4.12) to Mtot — 2.0 Mq. The total fermionic 
number is N ~ 920, which gives N ~ 4.7 x 10 57 through eq. (4.13). 

The potential U (</>) does not play any role in the solution we presented. The reason is that 
its value is much smaller than the pressure or the field derivative energy. The field essentially 
"rolls" on the effective potential U e tf = p. Only for f — > oo we expect U to become important. 
In this limit, our solution must be replaced by the time dependent cosmological solution. 

In fig. 2 we present the form of the solution if the chemical potential is [io/o~ ~ 0.95 while 
the remaining parameters remain the same. The deviations of 4>, A and B from 1 are much 
smaller than in the previous case. The fermionic gas is much more dilute, as its pressure is 
approximately 20 times smaller. The resulting astrophysical object has R ~ 22, R s ~ 0.41, 
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Figure 2: The solution of eqs. (4-14); (4-15), (4-16) for a model with a/GeV = 1, po/o~ ~ 0.95, 
4>i/M = 1, c = 1, C ~ 3.4 x 10 -47 . We plot (J), p, A, B, as a function of f. The object has 
R ~ 22, R, ~ 0.41, N ~ 10.5. 



which correspond to R ~ 7.4 km, Mtot — O.O47M0. The total fermionic number is N ~ 10.5, 
which gives JV~ 5.4 x 10 55 . 

The variation of the chemical potential results in a whole class of solutions depicted by the 
dashed line in fig. 3. We display the mass to radius relation Mtot{R)- This function has a 
maximum at R ~ 25 km, which separates two branches of the curve. The left branch is denoted 
by a thinner dotted line in fig. (3) because it corresponds to unstable configurations. For a 
given value of the total fermionic number N there are two solutions with different values of M tot 
[18, 19]. The dotted line corresponds to the branch with larger mass. The solutions of figs. 1 
and 2 are denoted by a circle and a square, respectively, on this plot. 

If the value of the scale a changes while <j>i is kept constant, the line Mt t{R) retains its 
shape, but is rescaled by an overall factor. This is obvious from the form of eqs. (4.14)-(4.19), 
in which a does not appear explicitly after the rescaling of eq. (4.8), as it has been incorporated 
in the various dimensionless parameters. If a solution with given values of a, (f>\ and no is 
known, another solution can be generated through the replacements a — > aa, <p — > (ft, r — > r /a 2 , 
p — ► a 4 p, p — > a 4 p and <j)\ — ► /Xq — ► a po, R — > R/a 2 , M to t — ► M tot /a 2 . In fig. 3 we present a 
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Figure 3: The function R S {R) for various models. 



class of solutions with 4>i/M=l, cr/GeV=3/2 (solid curve) that demonstrates this point. Each 
point of this curve can be obtained from a point of the curve with <pi/M=l, cr/GeV=l (dashed 
curve) by dividing the coordinates by a factor 9/4. Lowering a results in compact objects whose 
maximum mass can be significantly larger than the solar mass, while their average density falls 

4 

~ a . 

Changing the asymptotic value §\ of the field, while keeping the scale a constant, modifies 
the shape of the curve M tot (R). In fig. 3 we plot M tot (R) for <p 1 /M=3/2, cr/GeV=l (short- 
dashed curve). Comparison with the case </>i/M=l, <r/GeV=l (dashed curve) shows that the 
maximum mass is significantly reduced, while the mass is a decreasing function of the radius 
for the whole stable branch. The form of the curve is very reminiscent of that of conventional 
neutron stars. The reason is that for large (f>± the fractional deviation in the interior of the 
solution from the asymptotic field value is not very large. The fermionic mass is reduced but 
does not approach zero. This is similar to the case of nuclear matter, for which the effective 
reduction of the mass of a free fermion (~ 1 GeV) is of the order of the nuclear binding energy 
(~ 15 MeV). In fig. 3 we also depict the case with (f>i/M=3/2, cr/GeV=2/3 (dot-dashed curve). 
This has the same asymptotic fermionic mass as the case (f>i/M=l, cr/GeV=l (dashed curve). 

It is interesting that the upper part of the stable branch of the cases with 4>\/M = 1 has 
the same form as for fermion soliton stars [18], or fermion Q-stars [19], or strange stars [27]. 
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For these cases the fermionic mass approaches zero in the interior and the binding energy is 
large, similarly to the fermion stars. The mass is an increasing function of the radius. However, 
the lower part of the stable branch differs from that of the fermion stars. The mass becomes a 
decreasing function of the radius, similarly to the conventional neutron stars. 

5 Summary and conclusions 

In this work we considered some of the implications for astrophysical configurations composed of 
dark matter of the possible interaction with a scalar field that is responsible for the dark energy 
of the Universe. We assumed that at cosmological scales the scalar field has an almost constant 
expectation value. The time dependence associated with the evolution of the dark energy density 
was neglected, as we assumed that the relevant time scale is very long. The interaction between 
dark matter and dark energy was modelled by assuming that the the mass of the dark matter 
particles depends on the expectation value of the scalar field. We developed a formalism based 
on the Thomas-Fermi approximation, which assumes that in a local frame at every point in space 
the particles are described by a Fermi-Dirac distribution with position-dependent temperature 
and chemical potential. We derived the Einstein equations and the equation of motion of the 
scalar field. 

In regions of high number density of dark matter particles, the scalar field is shifted from its 
cosmological value. The shift is in the direction that reduces the dark matter mass, so that the 
total energy of the configuration is minimized. This deviation from the asymptotic value is the 
manifestation of the attractive interaction between the dark matter particles mediated by the 
scalar field. 

We discussed the implications for the dark matter halos that comprise the bulk of the matter 
outside the core of galaxies, up to distances of 100-200 kpc. Our approach is based on the 
assumption that the dark matter can be treated as a gas of thermalized particles. It may seem 
unlikely that the formalism we developed could apply to the dark matter halos, as the particles 
in them interact very weakly. However, we showed that our formalism is equivalent to the model 
of the isothermal sphere, which gives a simple explanation for the approximately flat rotation 
curves in the absence of the scalar interaction. The effective temperature in our approach is 
proportional to the dispersion (v 2 ,) of the velocity of the dark matter particles. 

The presence of the new attractive force does not modify the distribution of dark matter so 
as to destroy the approximately flat rotation curves. The main new effect is that the velocity of 
a massive object orbiting the galaxy outside its core is not of the order of the typical velocity 
of the dark matter particles, as in the conventional picture. Instead, it is reduced by a factor 
(1 + k 2 ) 1 / 2 , where k 2 = AM 2 [{dm® / 'dtp) / mo] 2 quantifies the dependence of the dark matter mass 
ttiq on the scalar field. If k 2 is large, the typical velocity of the dark matter particles can be 
significantly larger than the rotation velocity. The latter quantity is directly measurable, and its 
value is used in order to deduce the velocity of dark matter particles for dark matter searches. 

For k 2 ^ 1 the typical velocity of dark matter in our neighborhood of the Milky Way exceeds 
the rotation velocity significantly. The flux of dark matter particles towards a terrestrial detector 
is larger roughly by a factor (1 + k 2 ) 1 / 2 relative to the k 2 = case. As a result, the counting 
rates are increased by the same factor. Existing bounds on dark matter properties from direct 
searches can be extended to include the case of non-zero k. The bound on the cross section for 
the interaction of dark matter with the material of the detector is strengthened by the factor 
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(1 + k 2 ) 1 / 2 . Corrections to this simple picture are also possible through the velocity dependence 
of various contributions to the counting rates, such as nuclear form factors. A separate study is 
required in order to determine the implications for the various experiments. 

Using the same formalism we also derived solutions that describe denser compact objects 
composed of dark matter. In the examples we discussed, the dark matter gas has zero tem- 
perature and the stability of the configurations is provided by the degeneracy pressure. The 
particle mass is determined by the scalar field through a Yukawa interaction. The compact ob- 
jects resemble neutron stars, but are composed of dark matter, while their density can deviate 
significantly from nuclear density. Similarly to the case of dark matter halos, the potential of 
the scalar field plays a minor role in the structure of these objects. 

Their mass to radius curve has a stable branch whose shape depends on the asymptotic value 
of the scalar field in Planck units <f>i/M. For (fti/M i> 1 the mass of the compact object is a 
decreasing function of the radius, and the curve resembles strongly the one for neutron stars. 
The reason is that for large <f>± the fractional deviation in the interior of the solution from the 
asymptotic field value is not very large. The fermionic mass is reduced but does not approach 
zero. This is similar to the case of nuclear matter. For 4>i/M ^ 1 the upper part of the stable 
branch has the same form as for fermion soliton stars [18], or fermion Q-stars [19], or strange 
stars [27]. For these cases the fermionic mass approaches zero in the interior and the binding 
energy is large. The mass is an increasing function of the radius. However, the lower part of the 
stable branch differs from that of fermion stars. The mass becomes a decreasing function of the 
radius, similarly to neutron stars. 

The mass to radius curve also depends on the scale a that determines the mass of the dark 
matter particles. We assumed that the particle mass is generated by the scalar field through a 
Yukawa term, with a Planck suppressed Yukawa coupling a/M. For (fti/M = 0(1) and a = 0(1) 
GeV, the most massive astrophysical object has M to t = 0(1) Mq and R = 0(10) km. Changing 
a — ► aa leads to the rescaling R — > R/a 2 , M tot — > M to t/ct 2 . Reducing the scale a results in very 
massive, but dilute astrophysical objects. For a = 0(10) keV, we have M to t = 0(1O 10 ) M and 
R = 0(1O 10 ) km. These objects are similar to the supermassive neutrino stars hypothesized in 
ref. [20]. For a = 0(1O -3 ) eV the radius approaches the horizon size. Conversely, increasing 
a gives rise to very dense configurations of smaller size. For a = 0(1) TeV and <j>\ = 0(1) M 
the mass of an unbound dark matter particle is 0(1) TeV. The most massive compact object 
composed of dark matter particles has Mtot = 0(1O -6 ) Mq and R = 0(1) cm. 

The configurations we described in this work correspond only to a subclass of the possible 
solutions for the system of a Thomas-Fermi gas interacting with a scalar field in the presence 
of gravity. We considered an interaction only through the fermionic mass. Also, the novel 
astrophysical objects we described were assumed to have zero temperature. This assumption 
is probably unrealistic, as it means that during gravitational collapse the dark matter particles 
lose most of their momentum, apart from the amount required in order to satisfy the exclusion 
principle. 

Modifying some of our assumptions leads to different solutions. For example, a non-zero 
value for the temperature of the dark matter gas is expected to result in more dilute objects of 
larger size. Our aim here has been to demonstrate that the presence of an interaction between 
the dark matter and the dark energy has many more observable consequences, other than the 
modification of the cosmological evolution. One interesting property of our solutions is that the 
potential of the scalar field does not play any significant role. The reason is that the resolution 
of the coincidence problem requires the characteristic scale of this potential to be of the order 
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of the Hubble scale. As a result, the potential is negligible at length scales smaller than the 
horizon. 

Some of the ingredients of the model we considered can be fixed by establishing the correct 
theory underlying the dark sector. The form of the dependence of the dark matter mass on the 
scalar field, the value of the parameter a that sets the scale for this mass, the type of the scalar 
potential are typical examples. On the other hand, there are parameters that can be determined 
only by studying the dynamical mechanism that leads to the formation of astrophysical objects 
composed of dark matter. One such quantity is the temperature of the fermionic gas. This could 
be either an effective temperature of a weakly-interacting gas, as in the case of galaxy halos, or 
the real temperature of a thermalized gas. 

Addressing the dynamical problem is beyond the scope of this work. The crucial question 
is whether the gas of dark matter particles can lose momentum fast enough for dense objects 
to develop. The dark matter is expected to have only weak interactions, while its coupling to 
the scalar field is suppressed by the Planck scale. Despite this, it has been demonstrated that 
massive astrophysical objects can form under gravitational collapse [20, 28], even in the absence 
of the attraction mediated by the scalar field. For light neutrinos it has been shown in ref. [20] 
that matter can be expelled in a series of bounces during the collapse, so that a condensed object 
is left behind. In ref. [28] the case of neutralinos with a mass of 100 GeV has been considered. 
Through numerical simulations it has been shown that compact objects as large as the solar 
system and with a mass of the order of the Earth mass start forming at redshifts z ~ 60 and 
survive until today. 

The presence of an additional scalar attractive force increases the likelihood of formation of 
compact objects. The details of their structure (distribution on the radius-mass plane, density 
profile), the probability of their formation and the resulting abundance can be obtained only 
through a dynamical calculation that goes beyond linearized gravity within a specific theory 
of the dark sector. An analytical approach could make use of the model of spherical collapse, 
which has already been used for the study of structure formation for coupled dark energy and 
dark matter [29]. 

The most promising possibility for the detection of the astrophysical objects we discussed is 
through gravitational microlensing. The sensitivity depends crucially on the ratio of the object 
radius to the Einstein radius 



for gravitational lensing. The zero-temperature solutions that we presented have a mass to 
radius curve approximately given by eqs. (4.11), (4.12). For a £ 10~ 3 eV they satisfy R < R E , 
so that the corresponding astrophysical objects are expected to produce a detectable lensing 
signal. This must be contrasted with the much more dilute dark matter subhalos in the absence 
of an attractive force mediated by a scalar field [28]. These have size larger than the Einstein 
radius and are more difficult to detect. 
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